Is it really possible to grow isotropic on-lattice 
diffusion-limited aggregates? 

S. G. Alvesf and S. C. Ferreira Jr4§ 

fDepartamento de Fisica, Universidade Federal de Minas Gerais, CP 702, 
30161-970, Belo Horizonte, MG, Brazil 

fDepartamento de Fi'sica, Universidade Federal Vigosa, 36571-000, Vigosa, 
MG, Brazil 

Abstract. In a recent paper (Bogoyavlenskiy V A 2002 J. Phys. A: Math. Gen. 
35 2533), an algorithm aiming to generate isotropic clusters of the on-lattice diffusion- 
limited aggregation (DLA) model was proposed. The procedure consists of aggregation 
probabilities proportional to the squared number of occupied sites (fc^). In the present 
work, we analyzed this algorithm using the noise reduced version of the DLA model 
and large scale simulations. In the noiseless limit, instead of isotropic patterns, a 45° 
(30°) rotation in the anisotropy directions of the clusters grown on square (triangular) 
lattices was observed. A generalized algorithm, in which the aggregation probability 
is proportional to , was proposed. The exponent v has a nonuniversal critical 
value Vc, for which the patterns generated in the noiseless limit exhibit the original 
(axial) anisotropy for v < Vc and the rotated one (diagonal) for v > Vc- The values 
Vc — 1.395 ± 0.005 and Vc = 0.82 ± 0.01 were found for square and triangular lattices, 
respectively. Moreover, large scale simulations show that there are a nontrivial relation 
between noise reduction and anisotropy direction. The case v = 2 (Bogoyavlenskiy's 
rule) is an example where the patterns exhibit the axial anisotropy for small and the 
diagonal one for large noise reduction. 



PACS numbers: 61.43.Hv,05.40.Fb,05.10.Ln,05.50.+q 
1. Introduction 

Although proposed more than two decades ago, the diffusion- limited aggregation (DLA) 
model [T] constitutes an actual and challenging theoretical problem. This simple 
nonequilibrium growth process generates fractal structures with nontrivial scaling, which 
was subject of recent studies |21 El IH El- In spite of its simplicity, the DLA model 
was related to several physical and biological applications, such as electrodeposition [HI, 
viscous fingering [71, bacterial colonies P|, and neurite formation In the DLA model, 
particles released at points distant from the cluster execute random walks of unitary 
length. If a walker finds a site neighboring the cluster, it irreversibly sticks to this site. 
Both, on- and off-lattice versions of the model were intensively investigated pT H ITT l IT^. 

§ To whom correspondence should be addressed (silviojr@ufv.br). 
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The clusters generated by the on-lattice DLA model are unavoidable very sensitive 
to the lattice anisotropy ^21 The anisotropy is also a determinant factor present 
in other on-lattice growth processes as, for example, the Eden model jEl ^1 A 
procedure widely used to expose cluster anisotropy is the noise reduction, which consists 
of a set of counters at each growth site that are increased by a unit every time these 
sites are visited. An empty site is occupied only after its selection for M times. In figure 
H clusters generated by the original DLA model on the square lattice with distinct M 
values are shown. One can see that the effects of the axial anisotropy are enhanced as 
the noise reduction parameter is increased. 
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Figure 1. DLA patterns grown on the square lattice with distinct noise reduction 
parameters: (a) M = 2° (without noise reduction), (b) M = 2^, (c) M = 2^. The 
regions dehmited by the boxes correspond to 200 x 200 lattice unities. The number of 
particles in the clusters varies from 1900 (M = 2^) to 2800 (M = 2"). 

Recently Bogoyavlenskiy proposed an aggregation rule for the DLA model aiming 
to remove the anisotropy of clusters grown on square lattices |I7j. We refer to the 
Bogoyavlenskiy' s algorithm for the DLA model as BDLA along the text. The rule is the 
following. When a walker reaches an empty site neighboring the cluster, it will stick to 
this site with a probability Pk given by 

P,:P2:Ps:---:Pk = l':2^:3'---k', (1) 

where k is the number of occupied neighbors of the walker. When the aggregation 
does not occur, the particle is excluded and a new walk begins. This idea was based 
on a previous work, in which was verified that an aggregation probability proportional 
to the squared mean density is a special rule which generates isotropic patterns in a 
quasicontinuum extension of the DLA model 18,. Deviations from this rule (P oc p'^ 
with a ^ 2) result in anisotropic patterns with preferential growth along the axial or 
diagonal directions if a < 2 or « > 2, respectively. Three procedures were used to test 
the algorithm validity: (i) noise reduction algorithms with M < 2^; (ii) evaluation of 
the angular distributions of the particle density; and (iii) large scale simulations (10^ 
particles). 

In the present work, we analyze BDLA clusters grown in square and triangular 
lattices using large noise reduction and large scale simulations. The computational 
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procedures and the noiseless limit results are discussed in section |21 A generalized 
version of this algorithm is presented and discussed in section El and some conclusions 
are drawn in sectional 

2. Noise reduction in the Bogoyavlenskiy's algorithm 




Figure 2. BDLA clusters grown on square lattices using distinct noise reduction 
parameters: (a) M = 2° (without noise reduction), (b) M = 2^, (c) M = 2^, and (d) 
M = 2'^. The regions delimited by the boxes correspond to 200 lattice unities and 
the number of particles varies from 3000 (M = 2°) to 7000 {M = 2'^). These noise 
reduction parameters are the same used in reference |17| . 

In the original model proposed by Witten and Sander P, the particles are released, 
one at a time, at points distant from the cluster and execute random walks on a square 
lattice. If a particle reaches an empty site neighboring the aggregatef, it irreversibly 
sticks to this site. The distance between the center of the lattice and the hypersphere 
in which the walkers are released can be chosen as few lattices units larger than Rq, the 
maximum distance from the center of the lattice of a particle lying on the aggregate. 
Moreover, if the particle is too far way from the cluster, it is excluded and a new one 
released at the launching hypersphere. We used the launching and killing radius equal 
to Ri = Rq + 5 and Rk = 100i?o, respectively. Since this algorithm is very inefficient, 
the analysis with large noise reduction becomes prohibitive even for small systems. 
Therefore, we used a standard method in which the particles are allowed to execute 
long steps if they are far way from the clusters fOl ^1 QHl 120] • 

For BDLA clusters grown on the square lattice, the aggregation probabilities are 
given by Pk = where k = 1,2, or 3 is the number of occupied neighbors. 

Patterns corresponding to distinct noise reduction parameters M are shown in figure 
121 These patterns suggest the absence of anisotropy, since the cluster with the largest 
noise reduction (M = 2'*) does not exhibit preferential growth directions. Indeed, this 
can be reinforced comparing figure |2td) with the corresponding one for the original 
DLA model shown in figure ^c). Only the capillary width of these patterns, which 
enlarges with M, was modified. However, the noise reduction parameters used to 

t At the beginning of the simulation, a single particle at the center of the lattice constitutes the cluster. 
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Figure 3. BDLA clusters grown on square lattices using large noise reduction. The 
noise reduction parameters are (a) M = 2^, (b) M = 2^, and (c) M = 2* and the 
regions delimited by the boxes correspond to 200 x 200 lattice unities. 

generate these patterns are relatively small. In figure El simulations with larger M 
values are shown. For M = 2^, one can clearly observe that the diagonal directions have 
preferential growth. This becomes more evident as larger M values are used. In the 
noiseless limit, a regular fourfold structure with a 45° rotation in relation to the axial 
directions was found. Therefore, grounded in the later results one can conclude that 
the Bogoyavlenskiy's algorithm does not remove the lattice anisotropy. Instead, it only 
causes a rotation of the anisotropy direction. 

The physical origin of this rotation can be understood through the following 
argument. In square lattices, the net effect of the anisotropy is to increase the growth 
probabilities in the axial directions. Using an aggregation probability increasing with 
the number of occupied neighbor sites, the growth of the sites with two or three 
occupied neighbors is favored. Due to the ramified morphology of DLA clusters, empty 
sites with three occupied neighbors are much more rare than those with one or two. 
Thus, the effective result of this rule is to increase the growth probability of sites 
with two occupied neighbors and, consequently, to increase the growth probabilities 
in the diagonal directions. If the growth probabilities of empty sites with two occupied 
neighbors become larger than those of the sites with one occupied neighbor, the rotation 
of the anisotropy direction emerges. 

Another important feature of the BDLA model is the dependence on the lattice 
structure of the anisotropy rotations, as one can observe in figure lU where clusters 
grown on triangular lattices are shown. The aggregation probabilities for BDLA are 
given by 

= where = 1, ■ ■ ■ , 5. (2) 

In figuresEl^a) and (b) we show patterns generated by the original DLA model for M = 2 
and M = 2^. In turn, figures EJ^c) and (d) illustrate the corresponding patterns of the 
BDLA model. Similarly to the results of the square lattice, one can see that the effects 
of the algorithm are the enlargement of the capillary width and a 30° rotation in the 
anisotropy directions. In the noiseless limit, the patterns are characterized by a regular 
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sixfold structure in analogy with the square lattice case. 




Figure 4. DLA((a) and (b)) and BDLA ((c) and (d)) clusters grown on triangular 
lattices. The noise reduction parameter M = 2 was used in (a) and (c), while M = 2* 
was used in (b) and (d). The boxes correspond to 200 x 200 lattice unities. 

In order to confirm the results predicted by the large noise reduction analysis, we 
extend the large scale simulation performed by Bogoyavlenskiy ^Tj. In figure large 
aggregates generated using the DLA and BDLA models are shown. In figures El^a) and 
(b), we reproduced the simulations shown in reference in which clusters containing 
10^ particles were generated for both rules. In figure El^a), where the original model was 
used, the effects of the anisotropy are neat but, they do not appear in figure El^b). Thus, 
these figures suggest the anisotropy removal. However, if we continue the simulation 
shown in figureEfb) up to the cluster reaches 10^ particles, the axial anisotropy emerges, 
as one can see in figure Efc). 




(a) (b) (c) 

Figure 5. (a) Clusters with 10^ particles generated by the DLA model and clusters 
generate by the BDLA model with (b) 10^ and (c) 10^ particles. These clusters were 
grown in lattices 2000 x 2000 in (a) and (b), and 7000 x 7000 in (c). 
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3. The generahzed algorithm 

The distinct geometries observed for the noiseless DLA clusters grown using or not the 
Bogoyavlenskiy's rule suggest that this procedure can be generalized. Therefore, we 
modified the algorithm by defining an aggregation probability proportional to a power 
of the number of occupied neighbors, i.e., 

Pk = ( — 1 where k = 1, ■ ■ ■ ,n. (3) 

Here, n is the lattice coordination number and u a generalization parameter. For u = 
we recover the original DLA model and z/ = 2 the BDLA. Therefore, one expects a 
critical value Uc G (0, 2) delimiting these regimes. This aggregation rule is similar to 
that used in a modified DLA model investigated by Batchelor and Henry [22] • In this 
model, which was studied in square lattices, the aggregation probabilities are given by 
Pk = r^~^, where r G (0, 1] is a parameter. The noiseless limit of this model reveals a 
rich variety of morphologies that obey a nontrivial relation with the r parameter. Notice 
that in this aggregation rule the number of occupied neighbors appears in the power, 
contrasting the present generalization, in which this number appears in the basis. 




Figure 6. DLA clusters grown using the generalized algorithm. The noise reduction 
parameters are equal to M = 2^^ and M = 2^^ for square ((a)-(c)) and triangular 
((d)-(f)) lattices, respectively. The generalization parameters are (a) ly — 1.10, (b) 
ly = 1.36, (c) 1/ = 1.40, (d) ly = 0.77, (e) ly = 0.80, and (f) ly = 0.83. The boxes 
correspond to 200 x 200 lattice unities. 

Figures a)- (c) show clusters generated in square lattices for three values of the 
parameter u. Only the limit of large noise reduction was considered in this figure. For 
small u values, the morphologies are ruled by the usual square lattice anisotropy. When 
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the u values are increased, the clusters exhibit small diagonal modes coexisting with 
the fourfold patterns (e.g., figure 1(a)). As the u exponent is increased, the diagonal 
modes become more neat (e.g., figure 1(b)), up to a critical value z/c, when the diagonal 
anisotropy suddenly becomes the rule (figure 1(c)). Simulations with M = 2^^ (not 
shown in figure Ej) show that the critical value of u is in the interval 1.39 < Uc < 1.40 
and, therefore, the best estimate which can be made about the exponent value is 
Uc = 1.395±0.005. The influence of the u parameter in the morphology of clusters grown 
on triangular lattices is illustrated in figures intd)-(f). Even using a noise reduction 
16 times larger than that used for figures infa)-(c) (M = 2^''), these patterns do not 
display the closely regular structures observed for the square lattices. These simulations 
provided a critical u parameter in the range 0.81 < < 0.83 implying that the best 
estimate of the critical parameter for triangular lattices is Uc = 0.82 ± 0.01. 

The relation between the cluster anisotropy and the parameters u and M is 
nontrivial, as indicated in the morphological diagram shown in figure [7| These clusters 
were generated with lattices of size 2000 x 2000 and several values of u and M. It is 
clear from this diagram that for small u values the axial anisotropy is dominant while 
the diagonal one rules the patterns for large u. However, the dependence on the noise 
reduction is more complex. For u < = 1.395 ± 0.05, the noise reduction acts like in 
the DLA model, i.e., the larger the noise reduction more evident the axial anisotropy. 
Similarly, for u above an upper threshold u* ^ 3.4 ± 0.2, the noise reduction always 
enhances the diagonal anisotropy. But, for intermediate u values, the noise reduction 
enhances the axial anisotropy for small M while the diagonal one becomes the rule for 
large M. This is the case for u = 2 which corresponds to the BDLA model. In particular, 
the largest noise reduction used in the reference M = 2^, corresponds to a point 
near to the region of transition. In this region, both axial and diagonal anisotropics 
have approximately the same weight, and the patterns are apparently not sensitive to 
the lattice anisotropy. 

The determination of the transition line requires a quantitative measure of the 
anisotropy. Due to the fourfold symmetry of the patterns grown on square lattices, the 
fourth circular harmonic fOl 1201 



where the sum is done over all particles of the cluster and 9n is the angular position 
measured in relation to the origin, is an useful quantity since (cos 4^) > for 
patterns with axial anisotropy, (cos46') < for patterns with diagonal anisotropy, 
and (cos 46') = for isotropic or eightfold| patterns. Also, we evaluated the Fourier 
transform of the mass distribution defined by ||24j 




(4) 



n=l 




(5) 



J Axial and diagonal anisotropy are present with the same weight. 



Is it really possible to grow isotropic on-lattice diffusion-limited aggregates? 8 




M = 20 M = 2i M = 22 M = 2^ M = 2^ M = 2^ 

Figure 7. Morphological diagram for the generalized BDLA model using square 
lattices with 2000 x 2000 sites. The corresponding and AI values are shown at 
the left of the rows and at the bottom of the columns, respectively. The numbers of 
particles for these aggregates are of the order of 10^. 



Again due to the symmetries of the system, we are only interested in the ratio 
Ra = C^/Cq. In figures lEfa) and (b), we show (cos46') and R4 as functions of the 
parameter u for distinct system sizes and noise reduction M = 2. As expected, (cos 4^) 
is positive for small u values and negative for large values, vanishing for some value 
u'. The accuracy of these data allows us to estimate u' = 3.0 ±0.1 which, in this 
accuracy limit, is independent of the system size. The curves R4 against u exhibit 
positive minima at z/ = z/' that, within the error bars, are approximately independent 
of the system size. For u distant from u' (|z/ — z/'| ^0), R4 and |(cos46')| are increasing 
functions of the system size for a fixed z/, demonstrating the axial {u < v') or diagonal 
(z/ > v') anisotropy enhancement. However, in inset of figure |Hfb), the zoom around 
the minima shows that R^^iy ~ u') slowly decreases with the system size. So, it is 
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Figure 8. Measures of the anisotropies as functions of v for distinct system sizes and 
noise reduction M = 2. These quantities were defined in the text. In all simulations, 
the averages were done over 50 or 200 independent samples for |i/ — > 0.5 or 
liy — v'\ < 0.5, respectively. Figure|HIc) shows R4 as function of the system size for the 
off-lattice DLA model (o) and the generalized model at = i^' (♦) with the respective 
logarithm regressions. The inset of (b) shows a zoom around the minima of the curves 
i?4 versus v. 
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important to compare the minima of these curves with the R4 values obtained for off- 
lattice DLA clusters. For this task, we simulated 100 off-lattice DLA clusters of linear 
sizes 240, 480, 960, 1400, 3000 and 5000 and the results are shown in figure lEfc). As one 
can see, R4 decreases logarithmically slow even for the isotropic off-lattice DLA clusters. 
In turn, one expects that i?4 ~ for homogeneous structures, but this is not the case for 
DLA clusters, justifying the very low convergence to the asymptotic value. Also, figure 
ISfc) shows the minima of the curves in figure|Hfb) as a function of the system size. Again, 
the data suggest a logarithmic decay, but it is slower than that obtained for off-lattice 
DLA-model. This very slow decay indicates that the cluster anisotropy is minimal for 
u = u', probably reflecting their eightfold structures. Notice that we cannot claim that 
the patterns could be isotropic since i?4 is also null for exactly eightfold morphologies. 
Finally, we repeated the previous procedure for other noise reduction parameters and 
found z/' = 3.4 ± 0.2, 2.35 ± 0.05, 1.92 ± 0.02, and 1.79 ± 0.01 for M = 1, 4, 16 and 32, 
respectively. 

Baker and Ball |23j studied the anisotropy imposed by the square lattice in a DLA 
model where, beside the sticking to the nearest neighbors of the aggregate, the particles 
can also stick to the next nearest neighbors with a probability p. The parameter p 
controls the direction of the anisotropy, i.e., for p = the patterns exhibit the axial 
anisotropy and for p = 1 they exhibit the diagonal one. In this work, through real- 
space renormalization group, they found two stable and one unstable fixed points in 
the flux diagram p — M. The stable points {p = 0, M = 00) and (p = 1,M = 00) 
represent the axial and diagonal anisotropies, respectively, while the unstable one leads 
to eightfold structures. These results are in analogy with our generalization, which 
also seems to have two fixed stable points (z/ = 0, M = 00) and {i^ = 00, M = 00). 
Moreover, the qualitative morphological diagram and the quantitative characterization 
of the anisotropy suggest two distinct regions, in which only one of the anisotropy 
directions is present. The accuracy of the present simulations does not allow us to 
define the exact line delimiting these regions. Therefore, renormalization group analysis 
will be necessary to determine the exact flux diagram. In figure El we show a schematic 
representation of the morphological diagram, in which the increase of the z/' error bars 
as M — > 1 is illustrated by the hachured region. 

4. Conclusions 

In this work we revisited a procedure used to simulate isotropic on-lattice DLA clusters. 
In this procedure, the aggregation probabilities are proportional to squared number of 
occupied neighbors, i.e., Pk oc k"^. Contrasting the claims of this work, we demonstrate 
that this procedure implies in rotations of the anisotropy directions. More precisely, 
we found a 45° and 30° rotations for clusters grown on square and triangular lattices, 
respectively. Also, we studied a generalization of this method, in which Pk oc k'^. 
The exponent u has a nonuniversal critical value z/^ for which the patterns exhibit the 
usual (rotated) anisotropy directions for u < Uc {i' > z/c) in the noiseless limit. The 
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Figure 9. Schematic representation of the morphological diagram. The region 
depicted in gray (white) is ruled by the axial (diagonal) anisotropy. The hachured 
region represents the uncertainty around the line delimiting the two regions of the 
morphological diagram. Typical morphologies corresponding to distinct regions of the 
morphological diagram are illustrated. 

values obtained for the square and triangular lattices were Uc — 1-395 ± 0.005 and 
Vc — 0.82 ± 0.01, respectively. 

The hypothesis of the anisotropy rotations was confirmed through large scale 
simulations of clusters grown using square lattices. In turn, we found a nontrivial 
relation between cluster anisotropy and the parameters related to the generalizations 
and to the noise reduction. To be more precise, for intermediate v values, small 
noise reduction enhances the axial anisotropy while large noise reduction enhances the 
diagonal one. Also, using a quantitative analysis of the anisotropy, we determined the 
crossover value v' between axial and diagonal anisotropics for distinct noise reduction. In 
particular, for the case without noise reduction, we found u' — 3.4 ±0.2. It is important 
to mention that all implementation details were varied as, for example, the random 
number generators and the use or not of the optimizations. In all cases, the same results 
were obtained. Thus, our final concluding remark is that the possibility of growing 
isotropic on-lattice DLA clusters has not been demonstrated yet. In the best scenario, 
one can mask the anisotropy effects including additional anisotropy growth directions. 
As a last comment, we would like to mention that, concerning the computer time, 
the simulations of the on-lattice generalized model at the point of minimal anisotropy 
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